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Multi-poly-Bernoulli numbers and related zeta functions 


Masanobu Kaneko*and Hirofumi Tsumura^ 


Abstract 

We construct and study a certain zeta function which interpolates multi-poly-Bernoulli 
numbers at non-positive integers and whose values at positive integers are linear combinations 
of multiple zeta values. This function can be regarded as the one to be paired up with 
the ^-function defined by Arakawa and the first-named author. We show that both are 
closely related to the multiple zeta functions. Further we define multi-indexed poly-Bernoulli 
numbers, and generalize the duality formulas for poly-Bernoulli numbers by introducing more 
general zeta functions. 


1. Introduction 

In this paper, we investigate the function defined by 


( 1 ) 


. .,kr;s) 


mJo 


and its generalizations, in connection with multi-poly-Bernoulli numbers and multiple zeta values 
(we shall give the precise definitions later in §2). This function can be viewed as a twin sibling 
of the function ^{ki,..., kr‘, s), 


( 2 ) 


^{ki, ...,kr-,s) 


1 r 

ns) Jo e* - 1 


which was introduced and studied in [3]. The present paper may constitute a natural continua¬ 
tion of the work [3j. 

To explain our results in some detail, we first give an overview of the necessary background. 
For an integer A: E Z, two types of poly-Bernoulli numbers {B^^} and defined as 

follows (see Kaneko [20] and Arakawa-Kaneko [3], also Arakawa-Ibukiyama-Kaneko 131): 


( 3 ) 

( 4 ) 


Lifc(l - e p _ T 3 {k)n 

1-e-* X « 

n=0 

Lifc(l - e~*) ^ (j(k)n. 
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where Lifc(z) is the polylogarithm function defined by 


( 5 ) 




E 

m=l 



(N<i). 


Since Lii( 2 ;) = — log(l — z), we see that Bn^ (resp. Cn^) coincides with the ordinary Bernoulli 
number Bn defined by 


te 


e* - 1 


=y Bj- I 

n! 


n=0 




resp. 


e* - 1 


nl 

n=0 


(k) (k) 

A number of formulas, including closed formulas of Bn ' and Cn ' in terms of the Stirling numbers 
of the second kind as well as the duality formulas 


( 6 ) 

( 7 ) 





that hold for k,n G Z>o, have been established (see [2D1 Theorems 1 and 2] and [HI §2]). We 
also mention that Brewbaker |9] gave a purely combinatorial interpretation of the number Bn ' 
of negative upper index as the number of ‘Lonesum-matrices’ with n rows and k columns. 

A multiple version of B^^ is defined in [U p. 202, Remarks (ii)] by 


( 8 ) 


= . 

^ ' n=0 


where 

(9) 


^^kx,...,kr{z) — 


l<rm<-<mr '"'1 '''"2 


is the multiple polylogarithm. Hamahata and Masubuchi Hills] investigated some properties of 
Bn , and gave several generalizations of the known results in the single-index case. Based 
on this research, Bayad and Hamahata |8] further studied these numbers. Furusho [121 p. 269] 
also refers to ll8|). 

More recently, Imatomi, Takeda and the first-named author [18] defined and studied another 
type of multi-poly-Bernoulli numbers given by 


( 10 ) 

( 11 ) 


LiA:i,...,A:^(l - e *) ^ ^ 

1-e-* ^ ” n! 

n=0 

Lifci,...,fc^(l - e~*) ^ ^(ki,...,kr)'^ 


for ki,... ,kr € Z. They proved several formulas for and Cn^'""'^''\ and further gave 

ik\ k ) 

an important relation between C'p_ 2 ' and the ‘finite multiple zeta value’, that is, 


( 12 ) 


E 

l<mi<-"<mr<p 


1 


m 


ki 

1 ■ ■ 



^(kl,...,k ,— Ijfcr —1) 


mod p 


2 









for any prime number p. 

The function ([5]) for ki,...,kr € Z>i can be analytically continued to an entire function 
of the complex variable s € c m Sections 3 and 4]). The particular case r = k = 1 gives 
^(l;s) = Sec'S + !)• Hence ^{ki,..., kr; s) can be regarded as a multi-indexed zeta function. 
It is shown in [Ij that the values at non-positive integers of ^{k; s) interpolate poly-Bernoulli 
numbers Cm\ 

(13) = 


for k G Z>i and m € Z>o. And also by investigating ^{ki, ... ,kr] s) and its values at positive 
integer arguments, one produces many relations among multiple zeta values defined by 


(14) 


Ci^li ■ ■ ■ 1 ^r) 


E 

l<mi <-"<mr 


1 

mi'l ■ ■ ■ 


(= Lih,...T(l)) 


for li,... ,lr G ^>1 with Ir > 2 ([U Corollary 11]). 

Recently, further properties of ^{ki,kr] s) and related results have been given by several 
authors (see, for example, Bayad-Hamahata laiz], Coppo-Candelpergher [10], Sasaki [28], and 
Young [3^1. 

In this paper, we conduct a basic study of the function ([T]) and relate it to the multi-poly- 
Bernoulli numbers Bn ’' as well as multiple zeta (or ‘zeta-star’) values. Note that the only 
difference in both definitions m and ([5]) is, up to sign, the arguments 1 — e* and 1 — e * of 
Lifci,,,,,A:r('2) in the integrands. One sees in the main body of the paper a remarkable contrast 
between ‘H-type’ poly-Bernoulli numbers and those of ‘C-type’, and between multiple zeta and 
zeta-star values. We further investigate the case of non-positive indices ki in connection with a 
yet more generalized ‘multi-indexed’ poly-Bernoulli number. 

The paper is organized as follows. In §2, we give the analytic continuation of piki,... ,kr',s) 
in the case of positive indices, and formulas for values at integer arguments (Theorems 12.31 
and 12.5p . In §3, we study relations between two functions r]{ki,... ,kr] s) and ^{ki,... ,kr] s) 
fProposition 13.2() . as well as relations with the single variable multiple zeta function (Definition 
[26jand Theorem [ST]). We turn in §4 to the study of rj(ki,... ,kr;s) in the negative index case and 
give a certain duality formula for Bm (Definition 031 and Theorems 14.41 and 14.711 . We 

carry forward the study of negative index case in §5 and define the ‘multi-indexed’ poly-Bernoulli 
numbers for {ki, ..., kr) G (mi,..., m^) G Z>q and d G {1,..., r} (Definition 

[5TD . which include ([8]) and (HOD as special cases. We prove the ‘multi-indexed’ duality formula 
for them in the case d = r (Theorem 15.4p . 


2. The function r]{ki, ..., s) for positive indices and its values 
at integers 


2.1 Analytic continuation and the values at non-positive integers 

We start with the definition in the case of positive indices. 

Definition 2.1. For positive integers ki,... ,kr G Z>i, let 


r]{ki,... ,kr;s) 


1 /•°°,,-i Lifc„...,fc,(l-e*) 

ns)J, l-e* 
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for s E C with Re(s) > 1 — r, where r(s) is the gamma function. When r = 1, we often denote 
r]{k;s) by r]k{s). 

The integral on the right-hand side converges absolutely in the domain Re(s) > 1 — r, as is 
seen from the following lemma. 

Lemma 2.2. (i) For ki,... ,kr E Z>i, the function — e*) is holomorphic for t E C 

with |Im(t)| < TT. 

(ii) For ki,... ,kr E Z>i and t E M> 0 j we have the estimates 

(15) Lifc,,„„fc,(l-e*) = 0(r) (t^O) 

and 

(16) Lifci,...,fc,(l-e‘) = (t^oo). 


Proof. As is well-known, we can regard the function Li^^^ ^^^(a;) as a single-valued holomorphic 
function in the simply connected domain C \ [l,oo), via the process of iterated integration 
starting with Lii(z) = dz/{l — z). Noting that 1 —e* E [l,oo) is equivalent to Im(t) = (2j-|-l)7r 
for some j E Z, we have the assertion (i). 

The estimate ()15p is clear from the definition of Li^^ because its Taylor series at 

z = 0 starts with the term z^■ ■ ■ r^’’. As for (I16p . we proceed by induction on the ‘weight’ 
ki + ■ ■ ■ + kr as follows by using the formula 


(17) 


d j . ( \ _ j ^l-‘^fci,...,fcr-i,A;r—1('2) i^r > 1) 

I T^Li,,...,,,_,(.) {kr = l), 


which is easy to derive and is the basis of the analytic continuation of Lifc^^,,,^fc^(z) mentioned 
above. If r = fei = 1, then we have Lii(l — e*) = —t and the desired estimate holds. Suppose 
the weight k is larger than 1 and the assertion holds for any weight less than k. kr > 1, then 
by (flTl) we have 


|Lifci,(1 ^ ~ 


f 

f 




du 


U 


1 - 


;Lifci,...,fc,-i(l - e^)(-e^)(iu 


< 


/ 


,Lifci,„„fc,-i(l - e^) 


- 1 


dv + 


i: 


e^ -1 


(u := 1 - e^) 


LiA:i,...,fc^-i(l — e'^) 


dv 


for small e > 0. The former integral is 0(1) because the integrand is continuous on [0, e]. On the 
other hand, by induction hypothesis, the integrand of the latter integral is O as 

V ^ oo. Therefore the latter integral is O as t —>■ oo. The case of A:,- = 1 is similarly 

proved also by using (El), and is omitted here. □ 

We now show that the function r]{ki,... ,kr;s) can be analytically continued to an entire 

(ki k ) 

function, and interpolates multi-poly-Bernoulli numbers Bm ’ ’ at non-positive integer argu¬ 
ments. 
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Theorem 2.3. For positive integers ki,...,kr G Z>i, the function rj{ki,... ,kr', s) can be 
analytically continued to an entire function on the whole complex plane. And the values of 
rjiki,... ,kr',s) at non-positive integers are given by 

( 18 ) r]{ki,...,kr;-m) = (m € Z>o). 

(k) 

In particular, for k G Z>i and m G Z>o. 


Proof. In order to prove this theorem, we adopt here the method of contour integral represen¬ 
tation (see, for example, [3Ul Theorem 4.2]). Let S be the standard contour, namely the path 
consisting of the positive real axis from the infinity to (sufficiently small) e (‘top side’), a counter 
clockwise circle around the origin of radius e, and the positive real axis from e to the infinity 
(‘bottom side’). Let 


H{ki,...,kr]s) = [ F 

Je 1 


-1 Lifci,...,fc,.(1 e ) 


dt 


= (e^^*" - 


— e'- 

- ,_, Lifc,,...,fe^(l-e^) 

1 - 


dt+ t' 

JCe 


,-i Lifci,...,fc,(l - e*) 
1 — e* 


dt. 


It follows from Lemma r2.2l that H{ki,... ,kr;s) is entire, because the integrand has no singularity 
on C and the contour integral is absolutely convergent for all s G C. Suppose Re(s) > 1 — r. 
The last integral tends to 0 as e —>■ 0. Hence 


r]{ki, ...,kr\s) 


J'g27r2s 


1 

- i)r(s) 


H{ki,...,kr-,s), 


which can be analytically continued to C, and is entire. In fact r]{ki,... ,kr',s) is holomorphic 
for Re(s) > 0, hence has no singularity at any positive integer. Set s = —m G Z<o- Then, by 

COD, 


rj(ki,... ,kr; —m) = —-—^— H{ki,... ,kr; —m) 


27ri 
(-l)™m! 
2'ni 


« OO 

JCe n=0 


i-ty 


n\ 


■dt = 


This completes the proof. □ 

Remark 2.4. Using the same method as above or the method used in [4], we can establish the 
analytic continuation of ^{ki,... ,kr',s) to an entire function, and see that 

(19) = (m€Z>o) 

for ki,..., kr G Z>i, which is a multiple version of C2D- 

2.2 Values at positive integers 

About the values at positive integer arguments, we prove formulas for both f^{ki,... ,kr] s) and 
r/(fci ,... ,kr]s), for general index (ki,..., kr). These formulas generalize [U Theorem 9 (i)], and 
have remarkable similarity in that one obtains the formula for rj{ki,... ,kr]s) just by replacing 
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multiple zeta values in the one for ... ,kr',s) with multiple ‘zeta-star’ values. Recall the 
multiple zeta-star value is a real number defined by 


( 20 ) 


C{1 


1 , . . . , tj- 


Ir) = 


E 


h 


l<mi<---<mr 


• nir 


for li,... ,lr G Z>i with Ir > 2. This was first studied (for general r) by Hoffman in [16j . 

To state our theorem, we further introduce some notation. For an index set k = (ki, ..., k,.) G 
put k+ = (fei,..., kr-i, kr + 1). The usual dual index of an admissible index (i.e. the one 
that the last entry is greater than one) k is denoted by k*. For j = (ji,... ,jr) G Z> 0 ) we write 
|j| = ji + ■ ■ ■ + jr and call it the weight of j, and (i(j) = r, the depth of j. For two such indices k 
and j of the same depth, we denote by k-|-j the index obtained by the component-wise addition, 
k + j = (fei + Ji,..., fcr- -|- jr), and by 6(k; j) the quantity given by 



Theorem 2.5. For any index set k = (ki, ..., kr) G and any m G Z>i, we have 


(21) ^{ki,... ,kr-,m) = ^ 6((k+)*;j)C((k+)*-hj) 

|j|=m-l,(i(j)=n 


and 


(22) rj{ki,... ,kr]m) = (- 1 )^ ^ ^ 6((k+)*; j) C*((k+)*-h j), 

|j|=m-l,d(j)=n 

where both sums are over all j G Z>g of weight m — I and depth n := d{kf_) (= |k| -|- 1 — d(k)). 
In particular, we have 


f,{ki,..., kr',1) = C((k+)*) (= C(k+), by the duality of multiple zeta values) 


and 

V{ku...,kr-,l) = {-l)^-\*{{k+)*)- 

In order to prove the theorem, we give certain multiple integral expressions of the functions 
f,{ki,..., kr] s) and rj{ki,... ,kr] s). 

Proposition 2.6. Notations being as above, write (k+)* = (/i,...,f„). Then we have, for 
Re(s) > 1 — r, 

X --...^- dxi ■ ■ ■ dXn. 

gXlH- \-X„ _ -y gX2-\ - \-Xn _ I gX„ _ ^ 

V{ki,...,kr]s) = ^ [ ■[ {Xi + --- + Xn)"~^x[^~^---X^;j-^ 

lli=l i {h) ■ i ISj Jo Jo 

I ga:2H- hXn gX„ 

X -^-;- • -;-;...dXl • • • dXn- 

gXlH- \-Xn _ I gX2-\ - \-Xn _ gX„ _ 


6 













Proof. First write the index (ki ,..., kr) as 

[ki, ..., kr) = ( 1, ■ , 1 ,6i + 1,..., 1, ■ , 1 , bh + 1), 

ai — l cif^ — 1 

with (uniquely determined) integers h > 1, ai > 1 {1 < i < h), bi > 1 {1 < i < h — 1), and 
bh P 0. Then, by performing the intermediate integrals of repeated dz/{l — z) in the standard 
iterated integral coming from (|17ll . we obtain the following iterated integral expression of the 
multiple polylogarithm Liki,...,kriz)- 


w. u.)=r^r . 

Jo Xh Jo Jo Xh Jo ah' 


1 - ^ dxh-i 

1-Xh J Xh-l 



bh 


r'-l dXh-l 


r'*-! dxh-i 

Jo ^h—1 


lo Xh—1 , 


feh-'l-l 


n dX3 


r® dx3 

Jo X 3 


Jo X3 Jo 


63 — 1 


to 

1 — Xl ' 

dxi 

1 - X 2 , 

1 Xl Jo 


ah-P 


log 


1 - Xh-2 
1 - Xh-l 


<^h-l 


dXh-2 

Xh-2 


log 


1 — \ ^ dX2 dX2 

1 - X3 


dX2 dX2 dX2 

X 2 Jo X 2 Jo X 2 


dxi r 

Jo Xl Jq 


b2 — l 

ai 


dxi (-log(l - x))“^ dx 


ai' 


fei-i 


Here, to ease notation, we used the same variable in the repetitions of integrals Jq dxfx, and 
we understand Xh = z ii bh = 0. The paths of integrations are in the domain C \ [l,oo), and 
the formula is valid for z G C \ [1, oo). We may check this formula by differentiating both sides 
repeatedly and using (HZD. Putting z = 1 — e * and 1 — e*, changing variables accordingly, and 
suitably labeling the variables, we obtain 


ft ^4i+-+6^ ft2 

(23) Lifc,,...,fc,(l-e-‘)= / / •••/ — 

Jo Jo Jo e t 


H- i-bf^ — 


- t-bf^_i+2 _ f 


fe/i —1 


1 {tbi^ -+ i -hftji-i) 


a.h 


1 


1 


ah'- 


-+ l 


Jbi + --- + b^_l _ 2 


-I-6j,_2+2 _ -f 


bh-1 — ^ 


and 

(24) 


1 (4i+fe2 + l “ 


as' 


c>tb^+b2 + l — 


0 * 61+62 — 1 


0 * 61+2 _ I 


1 (4i+i —4i)” 
02 ! — 1 


0*61 _ I 


0*2 _ 


b2 — l 

1 1 

1 '^e^i - 1 ^^^1+-+^'*’ 


ft ^* 61 +...+ 6 ,, f^2 p- 

Li.i,...,..(l-e^) = (-l)^/ / ■■■/a 

Jo Jo Jo e “1 


bi-i 


**2 e *' Ji+-"+'>6 


iH- bbf^ — 


g*6iH-66^_i+2 

o*6lH-l-6^_i+2 _ 


bh-i 
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X 


QH+-- + bh-l 

g*S>lH- 'rb}^_l — 


- ybf^_2+2 

g^bjH-h6/,_2+2 _ 2 


1 fa+...+fc,_i+i - ■■■+'-'^- 1 +^ 

a/j! Qb^A -hb/,_l+l _ 2 


bh-^ —1 


1 (4i+b2+i-4i+62)“^eS+*>2+i eS+'>2 


o*fJl +2 


03 ! 


o^bj+b 2 + l — 


otbi+b2 — 


5*bl+2 _ I 


X 


1 eS 


62 —1 
e *2 1 


02 ! 


Jb^ + l _ 1 


Jb, _ I 


e *2 — 1 ai! e*i — 1 


dti clt2 ■ ■ ■ dtl)^-\ - \-bh‘ 


61-1 


The factor (—1)^ on the right of (IM |1 comes from (—= (—1)'’. Plugging (1^51) and 
into the definitions ([ 2 ]) and ([T]) respectively and making the change of variables 

t — Xi T ■ ■ ■ T Xfi, - \-bfi — X2 T ■ ■ ■ T Xfn - \-bfi — l — X^ + ■ ■ ■ Tx^, . . . , t2 — Xfi—l T Xni tl — Xn j 

we obtain the proposition. One should note that the dual index (k+)* = {li,... ,ln) is given by 

(k+)* = ( 1, • , 1 , ah + 1, 1, ■ , 1 , tth-i + 1,..., 1, ■ , 1 , ai + 1) 

bh bh-i — l fei —1 

and the depth n is equal to 61 + ■ ■ ■ + 6 /i + 1 , and that (the trivial) = 1 when k = 1. □ 

Proof of Theorem 12.51 Set s = m in the integral expressions in the proposition, and expand 
(xi + ■ ■ ■ + Xk)^~^ by the multinomial theorem. Then the formula in the theorem follows from 


the lemma below. 

Lemma 2.7. For li,... ,lr G Z>i with Ij. > 2, we have 

—1 —1 


□ 


1 fOO poo n-i- . . . r^Lr-l 

"WTW)L '"I 


— i e 


X2-\ - \-Xr _ \ 


gXr _ I 


dxi ■■ ■ dXr 


and 


Cih,---,lr) = 


00 ^.h-l , , . j.lr-1 


0 X 2 A - ^Xr 


nunij) 


10 JO 


gXlH- ^Xr _ I gX2^ - ^Xr _ I 


eXr — I 


dxi■■ ■ dXr 


Proof. The first formula is given in [H Theorem 3 (i)]. As for the second, we may proceed 
similarly by using n~^ = r(s)“^ dt to have 

00 CO ^ 

C{h,...,lr)= Y. E 

mi = l m 2 ,."?^r =0 


m''f (mi + m2)^2. (^mi + • • • + mrY^ 


1 /*oo roo 

Tr-rTnE E ■ 


-Ip-mixi ^ ^h-1 ^-{mi+m2)x2 , , , 


■ ■ - X. 


T-lg-(mi+-+m.K dxi---dXr 


fr~m S S 

^^3 — ^ ^ ^ ^ mi=l m 2 ,...,mr= 0 ' 


roo POO 


POO 

I —1 

I vL ^ cL 

Jo 


-1 . . . ™b-lp-mi(3:iH- ^Xr) . p-m2(x2^ -ha;,-) 


































nL, r(ii) 


r-f 




Jr-l 


• • • e 

- \-Xr 


dx^ 


■ dXr 


gXlH- VXr _ \ qX 2^ - \-Xr _ \ 


- 1 


-dxi 


dXr- 


□ 


We record here one corollary to the theorem in the case of r]k{'m) (compare with the similar 
formula in [U Theorem 9 (i)]). Noting (A: + 1)* = (1,..., 1, 2), we have 

fc-i 

Corollary 2.8. For k,m> 1, we have 

(25) Vk{m) = (jfc - l)C*(ji, • • • jfc)- 

JlH-= 


3. Relations among the functions rj and (^, and their con¬ 
sequences to multiple zeta values and multi-poly-Bernoulli 
numbers 


In this section, we first deduce that each of the functions r] and ^ can be written as a linear 
combination of the other by the same formula. This is a consequence of the so-called Landen- 
type connection formula for the multiple polylogarithm We then establish a formula 

for ^{ki,... ,kr',s) in terms of the single-variable multiple zeta function 


(26) 


C(^l) ■ ■ ■ ,lr\ s) — 


E 


h 




■ ml"' 


defined for positive integers li,... ,lj., the analytic continuation of which has been given in [3] 
(the analytic continuation of a more general multi-variable multiple zeta function is established 
in ID)- This answers to the question posed in §5 of [3]. As a result, the function rj{ki,... ,kr]s) 
can also be written by the multiple zeta functions of the type above. We then give a formula 
for values at positive integers of ^{ki,... ,kr] s), and hence of r]{ki,... ,kr] s), in terms of the 
‘shuffle regularized values’ of multiple zeta values, and thereby derive some consequences on the 
values of %(s). 

Let k = {ki ,..., kr) € be an index set. Recall that k is said to be admissible if the 
last entry kr is greater than 1, the weight of k is the sum ki + ■ ■ ■ + kr, and the depth is 
the length r of the index. For two indices k and k' of the same weight, we say k' refines k, 
denoted k ■< k', if k is obtained from k' by replacing some commas by -|-’s. For example, 
(5) = (2 -|- 3) ^ (2,3), (2,3) = (1 -|- 1, 2 -|- 1) ^ (1,1,2,1), etc. The standard expression of a 
multiple zeta-star value as a sum of multiple zeta values is written as 


C*(k)= C(k'), 

admissible 


where the sum on the right runs over the admissible indices k' such that k refines kb 

The following formula is known as the Landen connection formula for the multiple polylog¬ 
arithm ( [261 Proposition 9]). 
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Lemma 3.1. For any index k of depth r 


we have 


(27) Lik(^) =(-!)’■ Lv( 2 ). 

k^k' 

We can prove this by induction on weight and by using (fT7|l . see [26]. 

By using this and noting z/{z — 1) = 1 — e^ (resp. 1 — e~^) if z = 1 — e~^ (resp. 1 — e*), we 
immediately obtain the following proposition. 

Proposition 3.2. Let k he any index set and r its depth. We have the relations 


(28) 

r](k;s) = (- 

_l)7-i 





k^k' 

and 




(29) 

C(k;s) = (- 

.1)^-1 

Y 




k^k' 

Corollary 3.3. 

Let k be a positive integer. 

Then 

we have 

(30) 

Vk{s) = 

E 

C(k;s) 


\ii:weight k 


and 

( 31 ) Ck{s)= 

\<.:weight k 

where the sums run over all indices of weight k. Here we have written ^fc(s) for ^(k; s). 

Proof. The index (fc) is of depth 1 and all indices of weight k (admissible or non-admissible) 
refine (fc). □ 

We mention here that, also by taking k = (k) in Lemma (3.II and setting z = 1 — e* or 1 — e“*, 
one immediately obtains a kind of sum formulas for multi-poly-Bernoulli numbers as follows 
(compare with similar formulas in [171 Theorem 3.1]). 

Corollary 3.4. For k > 1 and m >0, we have 

(32) bW = (-1)™ 

^1 H- \-kr = k 

kj_ ,r>l 

and 

(33) CW = (_!)- ^ 

fclH- \-kr=k 

Next, we prove an Euler-type connection formula for the multiple polylogarithm. If an index 
k is of weight |k|, we also say the multiple zeta value C(k) is of weight |k|. 
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Lemma 3.5. Let k be any index. Then we have 

(34) Lik(l-2;)= ^ Ck(k'; - 2;)Lik'(2;), 

k',i>o ' 

where the sum on the right runs over indices k' and integers j > 0 that satisfy |k'| +j < |k|, and 
Ck(k';j) is a Q-linear combination of multiple zeta values of weight |k| — |k'| — j. We understand 
L\(ij{z) = 1 and |0| = 0 for the empty index 0, and the constant 1 is regarded as a multiple zeta 
value of weight 0. 

Proof. We proceed by induction on the weight of k. When k = (1), the trivial identity Lii(l — 
z) = Lii(l — z) is the one asserted. Suppose the weight |k| of k is greater than 1 and assume 
the statement holds for any index of weight less than |k|. For k = (ki,... ,kr), set k_ = 
{ki,..., kr-i, kr — 1) and k_|_ = (/ci,..., kr-i, kr + 1). 

First assume that k is admissible. Then, by dm) and induction hypothesis, we have 


d .r. . Lik_(l - z) 

—Lik 1 -z) = -^- 

dz 1 — z 


^ Ck_(l;j)Lii,...,i(l - ^)Lii(^ 


i,i>o 


the right-hand side being of a desired form. Here, again by (1171) . we see that 


-^Lii i(l - z)Lii(z) = ^ 
I — z dz 


^Lii,.,,,i(l - z)hYi+i{z 


i=0 


We therefore conclude 


Lik(l-2;) = - ^ Ck_(l; j) y~]Li^.^(l - z)ljYi+i{z) -b C 

i,i>o 


i=0 


with some constant C. Since lim^-^o Lii,,,., i(l — z)\j\\^ij^i{z) = 0, we find C = C(k) by setting 

j-i 

z ^ 0, and obtain the desired expression for Lik(l — z). 

When k is not necessarily admissible, write k = (ko,!,...,!) with an admissible ko and 


q > 0. We prove the identity by induction on q. The case g = 0 (k is admissible) is already 
done. Suppose ^ > 1 and assume the claim is true for smaller q. Then by assumption we have 
the expression 

Liko,i,...,i(l - 2 ) = ^ Ck'(m;j)Lii,...,i(l - 2;)Lim(2:), 

m,i>o 

where we have put k' = (ko, !,...,!). We multiply Lii(l — z) on both sides. Then, by the shuffle 

q-l 

product, the left-hand side becomes the sum of the form 

gLik(l-z)-F Lik^^i,...,i(l - z), 


kQiadmissible 


q-l 


and each term in the sum is written in the claimed form by induction hypothesis. On the other 
hand, the right-hand side becomes also of the form desired because 

Lii(l - 2 ;)Lii,.,.,i(l - z) = {j + l)Lii,.,,,i(l - z). 


j+i 


Hence Lik(l — z) is of the form as claimed. 


□ 
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With the lemma, we are now able to establish the following (see [H §5, Problem (i)]). 

Theorem 3.6. Let k be any index set. The function ^(k; s) can be written in terms of multiple 
zeta functions as 

?(k;s)= ^ Ck(k';j)r'^ ^V(k';s + j)- 

k',i>o \ J / 

Here, the sum is over indices k' and integers j > 0 satisfying |k'| + j < |k|, and Ck(k';j) is a 
Q-linear combination of multiple zeta values of weight |k| — |k'| — j. The index k' may be 0 and 
for this we set C(0; s + j) = C(s + j)- 

Proof. By setting z = e~^ in the lemma and using 


(35) 


Lii,...,i(^;) = 


(-log(l - z)y 


j! 


we have 


tP 


-t\ 


Lik(l-e *) = X] Ck(k'; j)—Lik'(e , 

k',j>o 

Substituting this into the definition ([2|) of ^(k;s) and using the formula ([H Proposition 2, (i)]) 


C(k;s) = f 

r(s) Jo 

we immediately obtain the theorem. 


oo ^s-l 


s* — 1 


Lik(e *) dt, 


□ 


Remark 3.7. This theorem generalizes [H Theorem 8], where the corresponding formula for 
Lii,.,.,i,fc(l - z) is 

r —1 


Lii,...,i,fc(l - z) = (-1) 


fc-i 


E 


Lii,...,i(l — z)Li 


ji, 


31H-|-3fe='r+fc 

V3i>l 




k-2 

r-l '~T^ 




As pointed out by Shu Oi, one can deduce Lemma [331 bv induction using [271 Prop. 5]. However, 
to describe the right-hand side of the lemma explicitly is a different problem and neither proof 
gives such a formula in general. See also [25] for a related topic. 


Example 3.8. Apart from the trivial case (!,...,!), examples of the identity in Lemma [3.51 up 
to weight 4 are: 


Li 2 (l - z) 
Li 3 (l - z) 
Lii, 2 (l - z) 
Li 2 ,i(l - z) 
Li 4 (l — z) 


-Li 2 ( 2 ;) - Lii(l - z)\Ai{z) + C(2), 

Lii,2(z) Li2,i(z) -h Lii(l - z)Ui^i{z) - C(2)Lii(2;) -b C(3), 
-Li3(2;) - Lii(l - z)U2{z) - Lii4(l - z)Ij\i{z) + C(3), 
2Li3(2:) -b Lii(l - z)U 2 iz) + C(2)Lii(l - z) - 2({3), 
-Lii,i,2(^) - Lii,2,i(2) - Li2,i,i(z) - Lii(l - z)lAi^i^i{z) 
+C(2)Lii,i(z)-C(3)Lii(z) + C(4), 
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Lii,3(l - z) 


Li 2 , 2 (l - z) 


Li 34 (l - z) 


Lii,1,2(1 - z) 


Lii,2,i(l - z) 


Lb,1,1(1 - z) 


Lii,3(5;) + Li 2 , 2 (^) + Lb, 1 ( 2 ;) + Lii(l - 2 ;)Lii, 2 ( 2 ;) + Lii(l - z)\A 2 ,i{z) 
+Lii,i(l - 2;)Lii,i(2:) - C(3)Lii(2:) + ^C(4), 

-Lb, 2 ( 2 :) - 2Li3,i(2;) - Lii(l - 2;)Li2,i(2;) - C(2)Lii(l - 2;)Lii(2;) 

—C(2)Li2(z) + 2C(3)Lii(2;) + ^C(4), 

-2Lii,3(2;) - Li2,2(2;) - Lii(l - 2;)Lii,2(2;) + C(2)Li2(2;) 

+C(3)Lii(l-z)-^C(4), 

-Li4(2;) - Lii(l - 2;)Li3(2;) - Lii,i(l - 2;)Li2(2;) 

-Lii,i,i(l - 2)Lii(2;) + C(4), 

3Li4(2;) + 2Lii(l - 2 ;)Li 3 ( 2 ;) + Lii,i(l - 2 ;)Li 2 ( 2 ;) + C(3)Lii(l - z) 
-3C(4), 

-3Li4(2;) - Lii(l - 2 ;)Li 3 ( 2 ;) + C(2)Lii,i(l - z) - 2C(3)Lii(l - z) 
+3C(4). 


Accordingly, we have 

i{2-s) = -C(2;s)-<(l;s + l) + C(2)C(s), 

e(3;s) = C(1,2 ;s)+C(2,1;s) + <(1,1;s + 1)-C(2)C(1;5)+C(3)C(s), 

^(l,2;s) = —C(3; s) — sC(2; s + 1)-—^C(l; s + 2) + C(3)C('S), 

e(2,l;s) = 2C(3;s) + <(2;s + l)+C(2K(s + l)-2C(3)C(s), 

^(4;s) = -C(l,l,2;s)-C(l,2,l;s)-C(2,l,l;s)-<(l,l,l;s + l) 
+C(2)C(1,1; S) - C(3)C(1; s) + C(4)C(s), 
e(l,3;s) = C(1,3 ;s)+C(2,2;s) + C(3,1;s) + <(1,2;s + 1) + <(2,1;s + 1) 

4——^C(i) 1; 'S + 2) - C(3)C(i; s) + -C(4)C(s), 

C(2,2;s) = -C(2,2;s)-2C(3,l;s)-<(2,l;s + l)-C(2)<(l;s + l) 
-C(2)C(2;s) + 2C(3)C(1;s) + ^C(4)C(^), 
e(3,l;s) = -2C(l,3;s)-C(2,2;s)-<(l,2;s + l) + C(2)C(2;s) 


+C(3)sC(s + 1) — ^C(4)C('S), 

■^(1)1; 2; s) = —C(4; s) — sC(3; s + 1)-—^C(2;s + 2) 

_ .(. + !)(.+ 2) 

b 

■^(1)2, l;s) = 3C(4; s) + 2sC(3; s + 1) H—^—^C(2; s + 2) + C(3)sC(s + 1) 
-3C(4)C(s), 

C(2,l,l;s) = -3C(4 ;s)-sC(3;s + 1)+C(2)^^^^C(5 + 2)-2C(3)sC(s + 1) 

+3C(4)C(s). 

From these and (l3Up of Corollary 13.31 we have for instance 

m{s) = ^(2;s)+C(i,i;s) 
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— “C(2; s) — sC(i; s +1) + C(2)C(s) h——^C( s + 2 ), 

mis) = e(3;s)+^(l,2;s)+^(2,l;s)+^(l,l,l;s) 

= c(3; s) + C(l) 2; s) + C(2,1; s) + sC(l) 1; 'S +1) ——^C(i; s + 2 ) 
“C(2)C(1; s) + C(2)sC(s + 1) H —^ -^C(s + 3), 

D 

mis) = e(4;s)+^(l,3;s)+^(2,2;s)+^(3,l;s)+^(l,l,2;s)+ai,2,l;s) 
+C(2,l,l;s)+e(l,l,l,l;s) 

= -C(4; s) - C(l, 3; s) - C(2, 2; s) - C(3,1; s) - C(l, 1,2; s) - C(l, 2,1; s) 
-C(2,1,1; s) - sC(l, 1,1; s + 1) + C(2)C(1,1; s) + 1; s + 2) 

-C(2)<(i;. +1) + C(2)^^C(. + 2) - £(i±il(i±^c(i;« + 3) 


6 


+Jc(4)C(.) + + 4). 


Before closing this section, we present a curious observation. Recall the formula 

^kim) = Cih_^,k + l) 

m—1 

discovered by Ohno [24]. Comparing this with the two formulas (12511 and [U Corollary 10], one 
may expect 

rjkim) = k + 1). 

m—1 

This is not true in fact. However, we found experimentally the identities 

(36) Vkim) = r]mik) 
and 

^ . . Ton - 

(37) 

J=1 

These are respectively analogous to the duality relation 

C(l ^k + l) = Cjl^^Ajm + l) 

m—1 k—1 

and the relation 

k-i oi-fc'i 


.y-i /2(l-2^ ^)C(/i;) (A:: even), 

= (^:odd). 


V(-i)^-k(i, ...,i,k-j + i) = r^^-"~^^^^^ 


which is a special case of the Le-Murakami relation [2^ (or one can derive this from the well- 
known generating series identity [2], m 

T{1- X)T{1-Y) 


1- c(i,...,i,A:-j + i)x"-^y^' 

k>j>l 


r(i - x-T) 


by setting Y = —X and using the reflection formula for the gamma function.) 

We are still not able to prove (ISGllPI . but could prove ([37)1 by using the following general 


^Quite recently, Shuji Yamamoto communicated to the authors that he found a proof. 
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formula for the value .^(k; m) and the relation 
‘height one’ multiple zeta values, see 


m 


Proposition 13.21 For other aspects of 


Proposition 3.9. Let k be any index and m > 1 an integer. Then we have 
(38) e(k;m) = (-!)— 


m—1 


where C™ stands for the ‘shuffle regularized’ value, which is the constant term of the shuffle 
regularized polynomial defined in m- 

Proof. By making the change of variable x = 1 — in the definition ([2]), we have 

C(k;s) = -^/' (-log(l-x))"“^Lik(x)—. 

^ [sj Jo ^ 

Put s = m and use (1351) to obtain 

dx 

^(k;m) = y Lii,.,,,i(x)Lik(x) —. 

m—1 


The regnlarization formula m Eq. (5.2)], together with the shuffle product of Lii^... ^ i(x)Lik(x), 

m—1 

immediately gives (f38]l . □ 

By using (1381) and (|30p . we can write ?/fc(m,) in terms of shuffle regularized values. The 
following expression seems to follow from that formula by taking the dual, but we have not yet 
worked it out in detail. 


Vkim) = 


? fm + k 


k 


C{m + k) — 


E 


h+- 


2<r<k+l 
■-{-jr=fTi + k — r — l 


H-h jr-1 

k — r + 1 


C(ji + 1, • ■ • ) jV-i + 1, jr + 2). 


4. The function r]{ki, ..., s) for non-positive indices 


In this section, as in the case of positive indices, we construct r/-functions with non-positive 
indices. It is known that \j\_k{z) can be expressed as 


Li-fc(^) 


P{z;k) 

(1 - z)fc+i 


for k € Z>o, where P{x] k) G Z[x] is a monic polynomial satisfying 


degP(x; k) 
x I P(x; k) 


1 (fe = 0) 

k {k > 1), 


(see, for example, Shimura [29l Equations (2.17), (4.2) and (4.6)]; Note that the above P{x;k) 
coincides with xPk+i{x) in [22]). We first extend this fact to multiple polylogarithms with 
non-positive indices as follows. 
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Lemma 4.1. For ki, ... ,kr € Z>o, there exists a polynomial P{x; ki,..., k^) G Z[x] such that 


(39) 

(40) 

(41) 


Li—( 2 ^) — 


P{z-,ki,. ..,kr) 

(]_ _ - hkr+r ’ 


degP{x;ki,... ,kr) = 
I P(x; ki,..., kr). 


ki-\ - 'rkr + r 


{ki = ■ ■ ■ = kr 
1 (otherwise), 


0 ) 


More explicitly, P{x] 0,0,... , 0) = 

'-V-' 

r 

Proof. We prove this lemma by the double induction on r > 1 and K = ki + ■ ■ ■ + kr > 0. The 
case r = 1 is as mentioned above. For r > 2, we assume the case of r — 1 holds and consider the 
case of r. When K = ki + ■ ■ ■ + kr = 0, namely ki = ■ ■ ■ = kr = 0, we have 


Lio,...,o(5;) = 


Z 


1 - z 


E 


E E 

mi<---<mi — 1 mr=mr—i + l 


yrUr-l _ . . . _ 




{i-zY^ 


which implies p^ - (l4T]l hold, and also P(x;0, ...,0) = x^. Hence we assume the case K = 
ki + ■ ■ ■ + kr — 1 holds and consider the case K = ki + ■ ■ ■ + kr{'> 1). We consider the two cases 
kr = 0 and kr > 1 separately. First we assume kr = 0. Then, by induction hypothesis, we have 




E 


fci k 

m-i ‘ ‘ ‘ 


-1 E 


— E 


1-z 




mr=mr—1 + 1 

kr —1 m,— 1 

V-1 


mi <-"<mr—1 


z P{z : ki,..., kr-i) 

— 2; - \-kr—l-\-T—l 


Let P{z; ki,..., kr-i,0) = zP{z; ki,..., kr-i). Then (f39]l - ()lTT) hold. 

Next we assume > 1. Then, using the same formula as in (|17p and the induction hypoth¬ 
esis, we have 


d f P{z : ki,... ,kr — Y'X 

^ dz V (1 ~ hfcr-l+r j 

_ z{P'{z : ki,... ,kr - 1)(1 - z) + {ki -h A:,. - 1 -|- r)P{z; ki,... ,kr - 1)} 

— 2;'jk,i-\ - \-kr-\-P 


If /ci = ■ ■ ■ = kr-i = 0 and kr = 1, then the numerator, that is, P(0,... , 0, —1) equals rz"^, using 
the above results. If not, the degree of the numerator equals /si A;,. -|- r — 1 by induction 

hypothesis. The both cases satisfy (j39p ~ (l4ip . This completes the proof of the lemma. □ 

Remark 4.2. In the case r > 2, P{x]ki,... ,kr) is not necessarily a monic polynomial. For 
example, we have Lio,-i(z) = 22;^/(l — zY, so P{x; 0,1) = 2x^. 
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We obtain from (1391) and (HOl) that 


(42) - e*) 


P{1 - e^;ki,... ,kr) 

g(fclH- \-kr+r)t 


0(1) = - = kr = 0) 

0(e“*) (otherwise) 


as t ^ oo, and from (HD) that 


(43) 




Therefore we can define the following. 
Definition 4.3. For ki,... ,kr G Z> 0 ) define 


(44) 


r]{-ki,..., 



hi—(1 e ) 


for s G C with Re(s) > 1 — r. In the case r = 1, denote rj{—k; s) by r]_k{s). 

We see that the integral on the right-hand side of (|44ll is absolutely convergent for Re(s) > 
1 — r. Hence r](—ki,..., —kr] s) is holomorphic for Re(s) > 1 — r. By the same method as in 
the proof of Theorem 12.31 for r]{ki,... ,kr] s), we can similarly obtain the following. 

Theorem 4.4. For ki,...,kr G Z>o, r]{—ki,...,—kr;s) can be analytically continued to an 
entire function on the whole complex plane, and satisfies 


rj{-ki,...,-kr;-m) = (m G Z>o). 


(45) 



It should be noted that f.{—ki ,..., —kr] s) cannot be defined by replacing {kj} by {—kj} in 
(j^. In fact, even if r = 1 and A: = 0 in ([2|), we see that 



which is not convergent for any s G C. Therefore we modify the definition ([2]) as follows. 
Definition 4.5. For ki,... ,kr G Z>o with (ki, ..., kr) (0,..., 0), define 



(46) 


for s G C with Re(s) > 1 — r. In the case r = 1, denote (,{—k] s) by ^_fc('S) for A; > 1. 

We see from ([42]) and (143]) that (]46]) is well-defined. Also it is noted that ^(A:i,..., A:,-; s) 
cannot be defined by replacing {—kj} by {kj} in (]45D for (kj) G 

In a way parallel to deriving Theorem 14.41 we can obtain the following. 

Theorem 4.6. For ki,... ,kr G Z>o with [ki,..., kr) (0,..., 0), ^(“^1) • • • > —^r] s) can be 
analytically continued to an entire function on the whole complex plane, and satisfies 


(47) 


e(-A:i,..., -kr] -m) = (m G Z>o). 


In particular, ^_fc(—m) = Cm {k G Z>i, m G Z>o). 
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Next we give certain duality formulas for Bn which is a generalization of ([6]). To state 

this, we define another type of multi-poly-Bernoulli numbers by 


r—1 


(48) 


a=0 


r — 1 


EMr 'T E 




a, / {h + ■ ■ ■ + Ir ~ 

T™-! . . . 

"I cO ^ 


I _ _ _ 


mi 


mi! • • • rrirl 


(k) (k) 

for s € C. In the case r = 1, we see that = Bm for fe G Z. Then we obtain the following 
result which is a kind of the duality formula. In fact, this coincides with Q in the case r = 1. 


v{-ki,...,-kr;s) = 


Theorem 4.7. For ki,... ,kr G Z>o, 

Therefore, for m G Z>o, 

(50) 


n(-fcl,...,-fcr) _ ai( 


Proof. We first prepare the following relation which will be proved in the next section (see 
Lemma 15.9h : 


(51) 


n 


- e*) 






fci ^, .^kf ^0 


holds around the origin. Let 


T{xi,...,Xr-,s) = ^ r]{-ki,...,-kr;s)-^^ 

fcl 5..., /Ct* ^ 0 

As a generalization of [181 Proposition 5], we have from (15111 that 

oo +S-1 _T_ _ gt) 


. . . .y.^r 

1 'T* 




2 ^OO 


1 - e^-=i^"(l - e*) 




El 1 — e“* (1 — e 


dt 




r(») 


a=0 

•oo 


mi j=l 


X 


r—1 


/ oo ' 

^s-lg(a-r)t 

i=i 


dt 


E(-i)‘ 


a=0 


r—1 


E 


nki 1-e 


=- Ei;=j 


a I ^^ (mi mr -|- r — aP 

mi,...,mr>0 


Therefore, by ()4S|1 . we obtain (14^ . Further, setting s = —m in ()4U]1 and using (I45|) . we obtain 

(pl) . □ 
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Remark 4.8. In the case r = 1, ()49ll implies r]_h{s) = Thus, using Theorem 14.41 we 

obtain the duality formula ([6]), which is also written as 

(52) r]-k{-m) = r]_rn{-k) 

for k,m ^ Z>o- This is exactly contrasted with the positive index case (I.S6I1 . Furthermore, 
by the same method, we can show that m) = C*!, ™ for A:,m G Z>o. Hence, using 

Theorem 14.61 in the case r = 1, we obtain the duality formula ([7]). 

Example 4.9. When r = 2, we can calculate directly from (H8|) that = 3“^ — 2“^. On the 

other hand, as mentioned in Lemma l4.ll we have Li_i_o(^) = ^^/(l “ Hence the left-hand 

side of (fTOll equals 

Li-i,o(l - e~*) ^ 1 -e~* ^ st _ 2 t 

l_e-t e-3i 

hence Bm = 3™ — 2"*. Thus we can verify Bm 


5. Multi-indexed poly-Bernoulli numbers and duality formulas 


In this section, we define multi-indexed poly-Bernoulli numbers (see Definition 15.111 and prove 
the duality formula for them, namely a multi-indexed version of ([6]) (see Theorem 15.411 . 

For this aim, we hrst recall multiple polylogarithms of *-type and of m-type in several 
variables defined by 


(53) 

(54) 




z 


mi 

1 


V 

l<mi<---<mr 


Si So 

1 "l2 


■ rUr 


E 

l<mi <---<mr 


E 


mi ^m2-mi 


mr—nir — l 
’ Z'p 




rrii m. 


Z 2 ' " 


rrir 

Ir 


—1 • • • (^1 IrY 

l,...,tr — 1 


for si,..., Sr G C and zi, • • • , Zr G C with \zj\ < 1 {I < j < r) (see, for example, [l3]). The 
symbols * and uj are derived from the harmonic product and the shuffle product in the theory 
of multiple zeta values. In fact, Arakawa and the first-named author defined the two types 
of multiple L-values L*{ki,..., kr; fi,..., fr) of *-type and {ki,..., kr‘, fi, ■ ■ ■, fr) of m-type 
associated to periodic functions {fj} (see [5]), defined by replacing {zy} by {fjim)} and setting 
(sj) = (kj) G on the right-hand sides of (f55|l and (fSTD for {ki,..., kr) G Note that 


(55) 


■ ■ ■ ’ ~ ^Li,...,Sr(n ■ ■ ■ ) Zr-lZr, Zr). 

i=l i=2 


Definition 5.1 (Multi-indexed poly-Bernoulli numbers). For si,..., Sr G C and d G {1, 2,... , r}, 
the multi-indexed poly-Bernoulli numbers are defined by 


(56) F{xi, . . . ,Xr]Si, ... ,Sr]d) = 


Li", 


l_e EU*-,...,i_e ^Dl-e- 


n^i 1- 
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^ oo 

= E 


nyi 1 

1 

1—1 


1 h (d) 

nyi 1 

(Ei=i c) 



E 


n{si,...,Sr),id.) 




mi,...,mr=0 


mi'. 


• mr 


where 6j{d) = 1 {j < d), = 0 {j > d). 

Remark 5.2. Note that ... ,z) = ^'iki,...,kriz) defined by ([9]). Suppose xi = ■ ■ ■ = 

Xr-i = 0 and {sj) = (kj) € in (l56]) . We immediately see that if d = 1 then 


jy{ki,...,kr),{l) _ T:){ki,...,kr) 


(see ®), and if d = r then 




_ Tri>(kl,...,kr) 


{m G Z>o) 


(m G Z>o) 


(see ([8])). 


Remark 5.3. Let 


(57) Ar = {(xi,..., Xr) G C*” I |1 — e ^'"| < 1 (1 < j < r)}. 


Then we can see that 


Li' 


Si ,Sr 


1 — e ^‘"=1 , 1 — e 


— Xr —1—a^r 


, 1 - e' 


(si ,... , §7- G C) 


is absolutely convergent for {xj) G Also F{xi ^..., Xr] 5i,..., 5^^; d) is absolutely convergent 
in the region x C^, so is holomorphic. Hence is an entire function, because 


r>(si,...,Sr),id) 




F{xi,...,Xr]Si,. 


J d) 




is holomorphic for all (si,..., G C'’. 

In the preceding section, we gave a certain duality formula for Bm " (see Theorem 14.711 . 
By the similar method, we can prove certain duality formulas for though they may 

be complicated. Hence, in the rest of this section, we will consider the case d = r. For emphasis, 
we denote by Bmi’/.’.’.’mr- Note that 6j{r) = 1 for any j. With this notation, we 

prove the following duality formulas. 


Theorem 5.4. For mi, ..., m^, ki,... ,kr G Z>o, 


(58) 


$(~kl,..., — kr) _ 


^kl,...,kr 


Now we aim to prove this theorem. First we generalize Lemma l4. II as follows. 

Lemma 5.5. For ki,... ,kr G Z>o, there exists a polynomial P{xi, ..., ki,..., kr) G Z[xi,..., x^] 
such that 


(59) 




nn;=i^.,n;= 


:2 ) 


, Zr—lZy, Zr,ki,...,kr) 


n 


j=i 


1 - n::= 


i/=j 
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(60) 


dega;jP(a;i,- • -.Xr-.ki, . . . ,kr) + I, 


V=J 


(61) (xi • • • Xr) I P{xi, . . . ,Xr]ki, . . . , kr). 

Set yj = n!)=j (1 < j < ^)- Then (l5^ implies 

P{yi,--- ,yr]ki,...,kr) 


(62) 


Li-fci,..._fc,(yi,---,yr) = 


nLi (1 - 2/j)^"^^ 


aE" fc.+i ■ 


Proof. In order to prove this lemma, we have only to use the same method as in Lemma 14.11 by 

induction on r. Since the case of r = 1 is proven, we consider the case of r > 2. Further, when 

iF = /ci + • • • + fcj. = 0, it is easy to have the assertion. Hence we think about a general case 

K = ki + ■ ■ ■ + kr{> 1). When kr = 0, we have 

Li-fci,...,-fc,(H, ...,Zr) = —_fc^_^(2:i, . . . , Zr-2, Zr-lZr) 

X Zj' 

Zr P(.Y\.j=l T • • • 1 —l^r j ) ■ ■ ■ ) l) 


1 — Zr 


nS(i-n::=,^ 


jp 


EEi k.+i 


Therefore, setting P{xi, ■ ■ ■ , ki, ..., kr-i,0) = XrP{xi, ..., x^-i; ki, ..., fer-i), we can verify 
(I59I)-(I6I]). 

Next we consider the case kr > 1. For k G Z>o, we inductively define a subset {c^m}o<j,i'<k+i 
of Z by 


(63) 

In fact, by 

and 


- f 

dz I ^ 


k-\-l k+1 


m^z^ = 


\m>l 


(1 — ^ 

^ ^ j=0u=0 




d_ 

dz 


E 

m>l 




m^z'^ = z 


dz 


\m>l 


-Lzm 


{k > 1), 


we can determine by (f63]l . Using this notation, we have 




d 


TPki,...,-kr(^^^ ■ ■ ■ ’ — Zr Li_^^^ _^^_|_]^(2;i, . . . , Zr) 


= Zr 


E 




mi • • • 1 ^ 


kr — l 7Tl\ VTlr 


mi<-"<mr—1 


1 


kj- kf 


(1 — Zr)’^^+^ 


EE 4 ! 


r-l 


(fer —l)^j + l 


^r-1 


E 


\^kr ^kr „(U —1) jy mr-l+j 
_l Z ^j=0 Z^u=0''j,u 'fl-r-l^r _ 

(1 — Zr)^’"*'^ 

ki fcr—1+^ mi / \mr—i 

m-i Z^ ^ ■ [Zr-lZr] u 


j=0 i/=0 mi<---<m ,— 1 

By the induction hypothesis in the case r — 1, this is equal to 


1 


(1 — Zr)’^^+^ 


kf kf 

EE 

j=0 v=0 


(kr-l) j + 1 


7=1 • • • ! Zr—lZr, ki, . . . , kr—2j kr—l T r*) 


n;=^(i-n::=,^ 


- Zr-lZr)>^r-l+u+l 
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Therefore we set 


j =0 i /=0 


, Xf—l^ /ci , . . . , kr—2 ? ^r— 1 “h ^) ■ 


Then this satisfies (f59]) ~ (l6T]l . This completes the proof. 

From this result, we can reach the following definition. 
Definition 5.6. For ki,... ,kr € Z> 0 ) define 


(64) 



□ 


for si,..., Sr S C with Re(sj) > 0 (1 < j < r). 

Lemma 15.51 ensures that the integral on the right-hand side of ()64p is absolutely convergent 
for Re(sj) > 0. By the same method as in the proof of Theorem 12.31 for r]{ki,... ,kr; s), we can 
similarly obtain the following. 

Theorem 5.7. For ki,... ,kr G Z> 0 ; v{~ki, • • •, —kr] si,..., Sr) can be analytically continued 
to an entire function on the whole complex space, and satisfies 

(65) hi-ki, • • •, -kr; -mi,-mr) = ■ ■ ■ ,mr e Z>o). 

Proof. As in the proof of Theorem 12.31 let 

(66) H{-ki,...,-kr;si,...,Sr) 



where S’” is the direct product of the contour C defined before. Note that the integrand on the sec¬ 
ond member has no singularity on C’’. It follows from Lemma fb.bl that H{—ki ,..., —kr', si,..., Sr) 
is absolutely convergent for any (sj) G C'’, namely is entire. Suppose Re(sj) > 0 for each j, all 
terms except for the first term on the third member of (1661) tends to 0 as e 0. Hence 

1 

'n{~kli • • • ) —kr', Si, . . . , Sr) = Yrr / o-wis ■ ivFT '\kP{ — kl, . . . , —kr'. Si, ..., s^), 

llj=i(e^ ^ - l)F(sj) 

which can be analytically continued to C^. Also, setting (si,..., s,.) = (—mi,..., —mr) G Z<q 
in (i66|) . we obtain (l6^ from ([M]) . This completes the proof. □ 
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Next we directly construct the generating function of rj{—ki, ..., —kr] si,..., Sr)- We prepare 
the following two lemmas which we consider when [xj) is in defined by (I57|) . 

Lemma 5.8. For (sj) G C'’ with Re(sj) > 0 (1 < j < r), 

(67) F{xi,... ,Xr-,si,... ,Sr]r) 


1 POO POO ^ ( 

'I 'i 








n ■ 

i=i 


— -•'{X — 

Proof. Substituting n“* = (l/r(s)) F~^e~'^^dt into the second member of (|56l) . we have 

"■ h-i 1 


F{{xj};{sj};r) = ^ 

— 1 j — 1 




n;=i ns, 


/ OO POO ' / -J \ ' 

• • • / n 1 (“£ ntj ] \ n 

j=l I \ U=1 / ) j=l 


We see that the integrand on the right-hand side can be rewritten as 


Ylf/ JJexp I I . 

i=i i=i 


u=, 


Hence we have 

F{{xj};{sj};r) = 




lj=l " \^3 

OO 

E n( 

l\ —1 j —1 

1 


n,^=ir(.,)yo 

1 


l_e-EE,-. 


-L nr' 

3=^ 


roo POO 


/o 


OO POO 


n< 

3=^ 

r 

n< 


Sj-l 

3 




1 - (1 - e 


-E:=,'^.)e-E:=,t 


■ n 

i=i 


gEEi 


_TT dt ■ 


n,^=ir(^.)yo 

This completes the proof. □ 

Lemma 5.9. Let zi^. .. ,Zr G C and assume that \zj\ (1 < j < r) are sufficiently small. Then 

( 68 ) 


n ^ — 


iL "I iL 


fcl,...,/L7'=0 




Set Zj = 1 — e^‘'=i (1 < j < r) /or (tj) G A^. Then 


(69) 


n 


gEi/=j ( X — 


e:=,' 4. 


i=i 


^ 1 — ( 1 — 6'^"=^ 






' 1 - e 

, . . . , -L C 


k\,...,kr=0 

In particular, the case t, = ■ ■ ■ = t^-i = 0 and tr = t implies (|5T]) . 


23 














Proof. We have only to prove ([68]) . Actually we have 


, , , rpkr 
1 y 




/C]^,... j/Ct* — 0 

oo oo r 

= E E n 

/ci,...,/Cr = 0 = l j = l 

oo r j 

= E n^rn' 

mi,...,7nr=l j=l fi=l 


hl---krl 

¥ 


E n 

j=l 


rr \ "O 


n 






j=i 1 - 


Thus we have the assertion. 

Using these lemmas, we obtain the following. 

Theorem 5.10. For ki,... ,kr € Z>o, 

(70) 77(-A:i,...,-fe,;si,...,s,)=«^"^’-’"^) 


Proof. By Lemmas 15.81 and 15.91 we have 
(71) F{xi,...,Xr;si,...,Sr]r) 


1 /•oo /•oo ' I 

rum)/ '"/ nj*? 


yj kf 




E 


oo /-OO 


'(1 - e 


e:;=u- 


■ JJ dtj 
i=i 


nt‘ 


10 JO 


nyir(sj; 

Li-A:i,...,-A:,(l - • • • > 1 - e*’') 


n-=i 




i=i 


ki\---kr\ 


□ 


for Re(sj) > 0 (1 < j < r). Combining (1^ . ([M|) and (ITT]) , we obtain ([70]) for Re(sj) > 0 
(1 < j < r), hence for all (sj) G C, because both sides of (f7n|) are entire functions (see Remark 

lOl) . □ 


Proof of Theorem \5.4[ Setting (si,..., s,.) = {—mi ,..., —m^) in (fTOjl . we obtain (f55|l from 
This completes the proof of Theorem 15.41 □ 

Example 5.11. We can easily see that 

T^Lii 1 .. ^ N ZlZ2{2 - Zl - Z 2 ) 

Li_i 0(^1 ?-^ 2 ) \2it ^ 0 ,— 1 (^1 ’'^ 2 ) 


{i-ziy{i-z2y 


{i-zini-z2r' 


Hence we have 


-l,o)^2-3^ = (2™ + 1)3" (m,n€Z>o). 


Therefore Bq = B^'^q = 3. Similarly we obtain, for example. 


4-1-2) 0) 


_ TOV _ -IQ toV ' 

1,0 ~ -^ 1,2 ~ -^ 1,2 


(-3-1) _ «{-!-2) 


3,1 ~ 1820, jlu2^2 


(-2-1) _ ini(-2-2) 


''2,1 


= 1958. 
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Remark 5.12. Hamahata and Masubuchi [141 Corollary 10] showed the special case of (|58l) 
namely 






(m,k G Z>o) 


(see Remark lOD . On the other hand, Theorem 14.71 corresponds to the case d = 1 r except 
for r = 1 (see Remark 15.21) . hence is located in the ontside of Theorem 15.41 Therefore, in (1501) . 
another type of multi-poly-Bernonlli numbers appear. 


Acknowledgments. The authors would like to thank Yasushi Komori for pointing out a 
mistake in an earlier version of this paper. 
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